On the accretion-induced QNM excitation of a Schwarzschild black hole 
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By combining the numerical solution of the nonlinear hydrodynamics equations with the solution of the linear 
inhomogeneous Zerilli-Moncrief and Regge-Wheeler equations we investigate the properties of the gravitational 
radiation emitted during the axisymmetric accretion of matter onto a Schwarzschild black hole. The matter 
models considered include quadrupolar dust shells and thick accretion disks, permitting us to simulate situations 
which may be encountered at the end stages of stellar gravitational collapse or binary neutron star merger. We 
focus on the interference pattern appearing in the energy spectra of the emitted gravitational waves and on the 
amount of excitation of the quasi-normal modes of the accreting black hole. We show that, quite generically in 
the presence of accretion, the black hole ringdown is not a simple superposition of quasi-normal modes, although 
the fundamental mode is usually present and often dominates the gravitational-wave signal. We interpret this 
as due to backscattering of waves off the non-exponentially decaying part of the black hole potential and to the 
finite spatial extension of the accreting matter. Our results suggest that the black hole QNM contributions to 
the full gravitational wave signal should be extremely small and possibly not detectable in generic astrophysical 
scenarios involving the accretion of extended distributions of matter. 

PACS numbers: 04.30.Db, 04.40.Dg, 95.30.Lz, 98.62.Mw 



I. INTRODUCTION 

It is well known within the framework of black-hole per- 
turbation theory H H H that quasi-normal modes (QNMs) 
(i.e., exponentially damped harmonic oscillations) dominate 
the gravitational-wave response of a non-spherically distorted 
black hole if the corresponding frequencies are part of the 
Fourier spectrum of the external source that moved the hole 
away from its equilibrium state. This perturbation is then ra- 
diated away in form of gravitational radiation until the black 
hole returns to its unperturbed, quiescent state. In practice, 
the QNMs excitation is triggered if the frequency of the per- 
turbing agent (e.g. an external matter source moving close 
to the black hole) is sufficiently close to the fundamental fre- 
quencies of the black hole, which then acts as an excited os- 
cillator. As a result, the QNMs and in particular the funda- 
mental mode (which is the one at the highest frequency and 
with the smallest damping time) represent the main feature of 
the gravitational waves emission only for a sufficiently com- 
pact perturbation; i.e., for sources whose characteristic scale 
is comparable with the width of the peak of the potential. 

These results are well-known since the early studies 
of Press J3l (see also Vishveshwara JH]), who considered 
the scattering of Gaussian gravitational-wave packets off a 
Schwarzschild black hole and noticed that the excitation of 
the QNMs of the black hole is more efficient for very narrow 
packets. Approximate relations to compute the efficiency of 
excitation of the various QNMs were introduced by Leaver [6] 
and by Andersson J3l for Gaussian pulses initial data with 
variable width. Since the QNMs spectrum is entirely de- 
termined by the black hole properties (i.e., mass, spin and 
charge), it is expected that the detection of the QNMs ring- 
down would provide a unique opportunity to unveil the phys- 



ical properties of a black hole. For this reason, the excitation 
of black hole QNMs have been studied in various astrophys- 
ical scenarios, such as the gravitational collapse of a (rotat- 
ing) neutron star to a black hole or the collision of two black 
holes. The presence of the QNMs in these situations has been 
confirmed either through the use of perturbation theory with 



im or through fully 



various degrees of sophistication 
relativistic numerical simulations 111 U 

However, under more general and realistic conditions, such 
as the excitation of the QNMs by accretion of matter, and 
that may be encountered in gravitational stellar collapse or bi- 
nary neutron star mergers, the gravitational-wave response of 
a black hole can be more c omp lex. For instance, the simula- 
tions performed by II 1611171 1 1811 showed that the gravitational- 
wave signal is not simply given by the superposition of 
exponentially-damped sinusoids at the QNMs frequencies and 
that, in some cases, the QNM ringing is only weakly excited 
and analysis in the frequency domain are needed 11911 . On 
the other hand, "backscattering" effects related to the slowly- 
decaying features of the scattering potential and interference 
effects turn out to play a crucial role for the correct interpre- 
tation of the results; a detailed discussion of these effects can 
be found in Ref. fll . 

In the case of a Schwarzschild black hole of mass M, we 
recall that the appearance of QNMs ringing is related to the 
peak of the curvature potential (also referred to as the Zerilli 
or Regge-Wheeler potential), that is located at r ~ 3 M; the 
backscattering, on the other hand, is related to the fact that the 
curvature potential decays as r~ 2 for r — ► oo. This behav- 
ior is responsible for the late time power-law tail £~( 2 ^+ 3 ) of 
the gravitational-wave signal, where I refers to the radiation 
multipole. 

Therefore, while early-time ringing is the result of a su- 
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perposition of exponentially damped sinusoids, and is domi- 
nated by the fundamental quasi-normal mode, at later times 
the ringing dies out and the signal is dominated by tail effects. 
However, in the transition from the "ringdown" to the "tail" 
phase, additional oscillations appear that cannot be attributed 
to any of the two regimes and that also seem to depend on the 
choice of the initial data (see 00). As we shall show in this 
paper for a broad sample of initial data, there could be inter- 
mediate regimes where the "ringdown" and the "tail" terms 
of the potential produce competing effects, so that the QNMs 
ringing and the backscattering effects can overlap, generat- 
ing complex waveforms. These effects were first noticed in 
Refs. lflq,[l7ll . although not discussed in detail there. 

In a recent paper lfl7ll . hereafter Paper I, a general analysis 
of the gravitational radiation emitted as a result of anisotropic 
accretion of matter shells onto nonrotating black holes and 
neutron stars was presented. That investigation made use of a 
procedure that combines the solution of the linearized Einstein 
equations for the metric perturbations with fully nonlinear hy- 
drodynamics simulations. Although the study of black hole 
perturbations produced by infalling matter has a long history 
and rich literature H2CX |2 lL 12211 . the approach outlined in Paper 
I proved to be useful for a number of reasons: i) it provided 
additional information on the black-hole's response to the dy- 
namics of point-like particles in the vicinity of black holes; 
ii) it helped understanding the basic black-hole's response to 
extended matter perturbations; Hi) it represented an effective 
way of studying black-hole physics in a linear regime without 
having to resort to full-scale numerical relativity simulations. 

One of the main results of Paper I was that, in the ideal- 
ized accretion processes considered, most of the energy is re- 
leased at frequencies lower than that of the fundamental quasi- 
normal mode (QNM) of the black hole, the spectrum consist- 
ing of a complex pattern, mostly produced during the accre- 
tion process rather than in the ringdown phase. More pre- 
cisely, the gravitational-wave emission was found to be dom- 
inated by a collection of interference "fringes" at frequencies 
of about a few hundred Hz, rather than by a single monochro- 
matic peak at the (higher) frequency of the fundamental mode 
of the black hole. Moreover, the width of these fringes was 
found to decrease rapidly with the initial position of the matter 
source. These results, which were already observed in other 
works 1 8, 23] in the case of a point-like particle falling onto a 
Schwarzschild black hole, also showed that the appearance of 
interference fringes in the energy spectra is much larger when 
the accreting matter is a shell of finite size and that the ef- 
ficiency in gravitational-wave emission is much reduced, be- 
coming almost two orders of magnitude smaller than in the 
case of point-like particles. An important feature of the cal- 
culations carried out in Paper I was the minimization of the 
initial gravitational-wave content; this turned out to be crucial 
to illustrate that the interference pattern was mainly due to the 
finite radial extension of the accreting source. 

The aim of the present paper is twofold. Firstly, we in- 
tend to complete the discussion started in Paper I on accreting 
quadrupolar shells onto a Schwarzschild black hole by extend- 
ing the parameter space of the initial models and by analyzing 
their impact on the gravitational-wave emission. In particu- 



lar we study how the energy emitted in gravitational waves, 
and the corresponding spectra, depend on the compactness of 
the shells as well as on their initial locations. In doing so we 
show that, for a finite-size source, the ringing of the black hole 
is much more complex than a simple superposition of QNMs 
and that the energy spectra (and in particular the interference 
fringes) are dependent on the choice of the initial data. 

Secondly, we improve the astrophysical relevance of our 
study by analyzing the gravitational radiation produced by 
thick accretion disks ll24l l25tl which accrete onto the black 
hole on dynamical timescales. We recall, in fact, that quasi- 
periodic oscillations of thick accretion disks (or tori) orbiting 
around Schwarzschild or Kerr black holes have been recently 
addressed as promising sources of gravitational waves 11811261 
IH HI IH in the kHz range. 

The paper is organized as follows: in Sec. [II] we review the 
theory of odd- and even-parity nonspherical perturbations of 
Schwarzschild spacetime, writing the inhomogeneous Zerilli- 
Moncrief and Regge- Wheeler equations in a form suitable for 
time-domain calculations. Sec. HTT1 briefly describes the nu- 
merical approach adopted for the simulations, while Sec. [IV] 
is devoted to the discussion of the results. Finally, Sec.lVlpro- 
vides a summary of the most important results and presents 
our conclusions. Unless otherwise specified, we choose ge- 
ometrized units (c = G = 1), and the black hole mass M is 
the unit of length. 

II. GAUGE-INVARIANT PERTURBATIONS OF THE 
SCHWARZSCHILD METRIC 

A. Odd and even-parity master equations 

The theory of gauge-invariant nonspherical metric pertur- 
bations of a Schwarzschild spacetime has a long history which 
has been recently reviewed in Refs. I3C& 13111 . We here simply 
recall that in this approach the spacetime metric is de- 
scribed by the background Schwarzschild metric g\ w plus a 
nonspherical perturbation h tliy . As a result of the spherical 
symmetry, can be expanded in three odd-parity and seven 
even-parity multipoles 

oo I 

9^=^+E E «") (0) + (C i ) (0) . (d 

1=0 m=-l 

with the odd multipoles transforming as (— and the 
even ones as (—1)' under a parity transformation (&, ip) — > 
(it — 9, 7r + ip). The presence of matter around the black hole 
is accounted through a "source" term in the linearized Ein- 
stein's equations, represented by a stress energy tensor i M „ 
which can also be expanded in multipoles. From the multi- 
poles of the perturbed metric, it is possible to build odd, M^L 

and even-parity, ^ < fj l , gauge invariant quantities which are so- 
lution of two equations, the Regge- Wheeler [32] and Zerilli- 
Moncrief IB3l 13411 equations, respectively and whose expres- 
sion in Schwarzschild coordinates is given by 

32^(0/0) „ 2 , x/o/cWo/c) _ „(o/c) m 
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Here the upper indices refer to the odd and even case re- 
spectively and = r + 2Mlog[r/(2A/) — 1)] is the "tor- 
toise" coordinate (see Refs. Oil 13511 for the most general 
covariant form of these equations). Explicit expressions for 



the scattering potentials and V^ w are well known and 
can be found, for example, in Ref. lT30h . The source terms 
in Schwarzschild coordinates can also be found in lfl7i l30ll . 
while the corresponding expressions valid for any coordinate 
slicing of Schwarzschild spacetime are reported in detail in 
Ref. JH. Once and are known, the gravitational- 
wave amplitude can be computed as 
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_2-f = -2 Y m (0, <P) are the 2 spin-weighted 
spherical harmonics [36] and the emitted power in gravita- 
tional waves is simply given by 
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Depending on the problem under investigation, we can ei- 
ther choose initial conditions that are time-symmetric (i.e., for 
which d t k£ m = d t Xim — 0) or not. In the first case, if the 
Hamiltonian constraint is satisfied, the momentum constraint 
is automatically satisfied too. In the second case, on the other 
hand, the momentum constraint must also be solved for dtki m 
if we provide a suitable ansatz for d r xtm- 

Whether it is possible to use time-symmetric initial data de- 
pends essentially on the astrophysical scenario under investi- 
gation and thus on the form of the stress-energy tensor of the 
matter, whose multipoles appear as sources in the Hamiltonian 
and momentum constraint equations. Clearly, time-symmetry 
is satisfied for matter configurations in equilibrium, or initially 
at rest and subsequently falling radially; in this case only the 
numerical solution of the Hamiltonian constraint is needed. 

This is the case, for example, of dust shells initially at rest at 
a finite radius. Because this configuration is intrinsically non 
time symmetric, it will invevitably produce spurious radiation. 
A way to minimize the latter is to freeze the sources of the 
perturbation equations up until the radiation has left the grid. 
This is the approach we will follow in most of our simulations. 



where the overdot refers to a derivative with respect to the 
Schwarzschild coordinate time. 



B. Initial data 

Specifying suitable initial data for nonspherical metric per- 
turbations taking place in astrophysical events is not trivial. 
The standard approach, in the case of even-parity perturba- 
tions, exploits the fact that can be written in terms of 
two gauge-invariant multipoles, namely the perturbed confor- 
mal factor k( m and the gravitational-wave degree of freedom 
Xim as introduced, for instance, in Eq. (8) of Paper I and 
whose relationship with the Zerilli-Moncrief equation is given 
by Eq. (18) of the same reference lf30L H3l As a result, the 
perturbed Hamiltonian equation and the equation for the mo- 
mentum constraint, i.e., Eqs. (45) and (46) of |3()|. provide a 
system of coupled ordinary differential equations for the un- 
knowns (ki m , Xim, dtki m , and dtXim)- The indetermination 
inherent in the solution of this system can be overcome, for 
example, by assuming xim = /3ke m B23I1 . where (3 > 0, so 
that the Hamiltonian constraint simply reads 
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where A = £(£ + 1) and Tq™ are the multipoles of the en- 
ergy density too of the matter source according to the mul- 
tipolar decomposition given by Eq. (13) of Paper I (see also 
Ref. 113711 ). The case f3 = corresponds to conformally flat 
initial data, since in this case the perturbed metric in isotropic 
coordinates is the Schwarzschild one modulo a conformal 
transformation. 



in. NUMERICAL FRAMEWORK 

The numerical approach adopted in our simulations is the 
same as the one described in Paper I. We only recall here that 
our hybrid approach implies the solution of both the nonlinear, 
relativistic hydrodynamics equations on a fixed Schwarzschild 
background, in the two spatial dimensions r and 9 (since we 
restrict ourselves to axisymmetric flows), and of the Regge- 
Wheeler and Zerilli-Moncrief, in the radial direction r*. Such 
"test-fluid approximation" is valid as long as the mass of the 
accretion flow is much smaller than the mass of the black hole; 
in our simulations this mass ratio is < 10% 

The hydrodynamics equations, once cast in conservation 
form, are solved using Godunov-type methods based on Rie- 
mann solvers (see e.g. l38ll and references therein), while 
Eqs. (ffj can be written as a first-order hyperbolic system over 
a one-dimensional (ID) grid expressed in terms of the radial 
coordinate r* and solved with some standard finite-difference 
method such as Lax-Wendroff (see Paper I for details). 

The computational domain used for the hydrodynamical 
evolution consists of N r x Ng grid-points, geometrically dis- 
tributed along r [that is, with grispacing Ar^+i = aAr^ with 
a > 1] and uniformly distributed along 9. When we evolve 
quadrupolar shells of dust plunging from large distances we 
set N r = 2000 and Ng = 20, while for the evolution of 
fluid accretion disks we use N r = 300 and Ng = 150 to 
reach the desired truncation error. The hydrodynamical grid 
extends from a point slightly outside the event horizon up to 
a finite radius, that depends on the initial data, much smaller 
than the extraction radius. At the inner boundary we impose 
inflowing boundary conditions, while at the outer one we fix 
the conditions of a tenuous stationary spherical atmosphere. 

The radial grid used for the time evolution of the Regge- 
Wheeler and Zerilli-Moncrief equations, on the other hand, is 
much more extended in radius than the hydrodynamical one 
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and overlaps with the latter. To avoid spurious boundary ef- 
fects, the ID grid typically extends from -2000 M to 5000 M 
and is covered by ~ 10 5 cells. At the inner and outer boundary 
of this grid we use standard Sommerfeld outgoing conditions. 
We stress that, due to the long-range action of the potential, 
these conditions are not completely non-reflecting, which is 
the reason why we still need large grids to properly capture 
the late time power-law decay of the waveforms (see below). 
We recall in passing that exactly non-reflecting boundary con- 
dition for the Zerilli-Moncrief and Regge- Wheeler equations 
have been recently proposed in ll39ll . 

While the gravitational waveforms are extracted using both 
the Regge- Wheeler and Zerilli-Moncrief equations, we also 
use as a comparison the Newtonian quadrupole formula to 
compute the amount of gravitational waves emitted. We use 
the modified definition of the quadrupole moment proposed in 
Ref. H 

and which reads 

I = 27T J dz J (hDW -p) Qz 2 -0 r 4 dr , (6) 

where D = pW is the conserved rest-mass density, W is the 
relativistic Lorentz factor, h is the relativistic enthalpy, p is 
the pressure, and z — cos 8. Expression © represents a slight 
modification of the Newtonian standard quadrupole formula 
(SQFi) and it was shown to provide the best agreement with 
the gauge-invariant waveforms when compared with alterna- 
tive definitions 11811 . 

IV. RESULTS 
A. Effects of the gravitational potential 

As mentioned in the Introduction, for a certain set of ini- 
tial conditions the ringdown phase of the gravitational-wave 
signal from a perturbed black hole shows the presence of 
the first QNMs of the black hole (typically, the fundamental 
mode and the first overtone) together with "backscattering" or 
"taif'effects. While the less damped modes are related to the 
interaction of the perturbations with the peak of the black hole 
Regge-Wheeler potential, the tail-effect cannot be associated 
to QNMs and are instead due to the long-range features of the 
potential and in particular to its slow decay with radius. 

In order to single out the contribution of the peak of 
the curvature potential on the generation of the resulting 
gravitational-wave spectrum, while suppressing the effect of 
the tail term, we solve a simplified version of Eq. (f2]i in 
which the Regge-Wheeler potential is replaced by the so- 
called Poschl-Teller potential i40ll4lll . which is exponentially 
decaying for r* — > ±oo. We therefore evolve in time an equa- 
tion of the type 

d^ e -d 2 r y e + u f y e = o, (7) 

where LTg is given by 

umax 

U e = 2 " ( g ) 

cosh [a(r* — r™ ax )] 



and r™ ax is the position of the peak of the Regge-Wheeler po- 
tential vy nax = Vg °\r™ ax ) and a is determined through the 
matching of the second derivatives at r™ ax ; i.e., d 2 Ui/dr 2 = 

d 2 V^ I dr 2 at r* = r™ ax . In the following, for simplicity we 
will consider only the I = 2 odd-parity perturbations. 

The Poschl-Teller potential is shown in the top left panel 
of Fig.Q]together with ■ Although the exponential decay 
of the former is evident for large r*, it is also apparent that it 
provides a good approximation to the Regge-Wheeler poten- 
tial near the peak. The accuracy of the approximation was first 
studied in Ref. B42I1 . where it was pointed out that the QNMs 
of the Poschl-Teller potential can be computed analytically. 
The frequencies found for the lower modes agree within a few 
percent with those obtained using the true curvature potential 
, computed numerically in 1I43I1 . 

We solve Eq. (O using a Lax-Wendroff method on an 
evenly spaced grid with Ar* = 0.1 M and with initial 
data given by a Gaussian pulse ^2 = N exp [—(r — r c ) 2 /a 2 ] 
where r c is the initial position of the pulse and N a normal- 
ization constant. The initial pulse is considered to be purely 
ingoing (i.e., with = d r ^2)- We compare the sig- 

nal extracted at r Q b s = 200 M for different values of a us- 
ing either the Regge-Wheeler or the Poschl-Teller potential. 
This is summarized in Fig.[T| which shows the waveforms ob- 
tained for different values of r c and a and expressed in the 
retarded time, u — t — r° hs . Note that the waveforms from 
the Poschl-Teller have been shifted in time so as to overlap the 
maxima and the minima of the ringdown phase of the Regge- 
Wheeler potential. This is necessary because the scattering of 
the pulse starts earlier for the Regge-Wheeler potential than 
for the Poschl-Teller one. 

Let us focus first on the top right panel of Fig. [T] which 
refers to a very narrow pulse (a = M) initially located at r c = 
100 M and shows with a solid and a dashed line the wave- 
forms obtained with the Regge-Wheeler and Poschl-Teller po- 
tentials, respectively. Clearly, the ringing is very similar in 
the two cases both in the wavelength and in the amplitude, 
but differences become apparent after u/M ~ 180, when the 
tail term of the Regge-Wheeler potential becomes dominant 
and the gravitational-wave signal is driven by the backscatter 
due to the r~ 2 decay and asymptotes the expected late-time 
decay according to Price's law [44]. Since the Poschl-Teller 
potential decays exponentially, no tail effects are found and 
the signal is still given by a superposition of exponentially 
damped harmonic oscillations. Smaller differences between 
the waveforms are however present also in the early part of the 
waveform (cf., the inset of the top right panel of Fig. [D and 
are, again, due to the fact that the interaction of the perturba- 
tion with the Regge-Wheeler potential starts "earlier" (i.e., at 
larger radii) than in the case of the Poschl-Teller potential. 

The two bottom panels of Fig. Q] refer instead to a wave- 
packet that is initially closer to the black hole (i.e., r c = 
50 M) and show the impact on the waveforms of an increas- 
ing width of the wave-packet. Most notably, the two panels 
show that as a is progressively increased, the effects due to the 
slow-decay of the Regge-Wheeler potential become progres- 
sively more pronounced, with the ringdown lasting progres- 
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Figure 1 : Comparison between the outcome of a scattering problem over the Regge-Wheeler (solid lines) and the Pbschl-Teller (dashed lines) 
potential. The waveforms are extracted at r b s = 200 M and refer to a choice of the parameters given by (from top to bottom and from left to 
right): a = M and r c — 100 M, a = 9.5 M and a = 11.5 M with r c = 50 M. The wider the pulse, the larger the differences between the 
two solutions. The top left panel shows a comparison between the profiles of the Regge-Wheeler and the Pbschl-Teller potentials. 



sively less and with the backscattering being correspondingly 
anticipated. For example, when a — 9.5 M (left-bottom panel 
of Fig.Q~|i and despite the fundamental mode is still recogniz- 
able, other oscillations are present already after u/M ~ 70, 
while for a = 11.5 M (right-bottom panel) the global wave- 
form is essentially overwhelmed by backscattering, the QNMs 
are almost absent and the signal is dominated by the "tail". An 
additional increase in the pulse width would make the QNMs 
disappear completely. Clearly, this behavior is not present in 
the case of a scattering off a Poschl-Teller potential, which 
is much less sensitive to the finite-size of the perturbation as 
long as it is smaller than the scale-height set by the exponen- 
tial decay. 

Although just a toy -problem, this simple comparison of the 



black-hole response when modeled with the Poschl-Teller po- 
tential is very useful to clarify that in the case of general ini- 
tial data with a finite size, the shape of the curvature potential 
strongly affects the gravitational waveform. In the next sec- 
tions we will refine this finding by considering more realistic 
perturbations such as those produced by the accretion of ex- 
tended distributions of matter. 

B. Quadrupolar shells of dust 

Following the analysis presented in Paper I, we here con- 
sider a number of aspects of gravitational-wave emission re- 
sulting from the accretion of quadrupolar shells of dust which 
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Figure 2: The 1 = 2 gravitational Zerilli-Moncrief function extracted at r Q b s = 500 M for a dust shell falling from ro = 15 M with 
a/M = 0.050. Left panel: the complete waveform. Right panel: investigation of the presence of the black hole QNMs: the dashed line was 
obtained through a non-linear fit with the first five modes and the comparison with the computed waveform is shown on both a linear and a 
logarithmic vertical scale. 



were not investigated in detail before. More specifically, we 
will here discuss: z) a more detailed analysis of the black 
hole ringdown phase, focusing on black hole QNMs and on 
backscattering effects related to the non-exponential decay of 
the black hole potential; it) the mechanism responsible for the 
production of the interference fringes in the energy spectra; 
Hi) the effect of using conformally-flat initial data; iv) a more 
detailed analysis of the energy released in gravitational waves. 
As in Paper I, the rest-mass density is parameterized as 

p = p + p exp[-(r -r ) 2 /cr 2 ]sin 2 z? , (9) 

where ro is the initial position of the "center" of the shell and 
(j controls its compactness. The background rest-mass den- 
sity p is chosen to be very small (i.e., ~ 10~ 22 ) to simulate 
the vacuum outside the black hole, while the normalization 
constant po is obtained from the condition that the volume in- 
tegral of Eq. (O gives p, the total mass of the shell, which we 
choose to be p = 0.01M. As we mentioned in Sec. Ill B I (and 
Paper I), the initial profile for ^P^o * s obtained after solving 
the Hamiltonian constraint for &2o with (3 = 0. Furthermore, 
to minimize the impact of surious radiation and produce initial 
data that is a "almost" time-symmetric, the shell is kept frozen 
at ro (i.e., the hidrodynamics equations are not evolved) up 
until the spurious initial gravitational-wave pulse leaves the 
numerical grid. We shall show at the end of this Section how 
waveforms and energy spectra change when "genuine" time- 
symmetric and conformally flat initial data are implemented. 

We start by considering the plunge from ro = 15M with 
a/M ~ 0.050 (this particular model was discussed also in 
Paper I) and whose complete waveform is shown in the left 
panel of Fig. [2] As expected, the gravitational-wave sig- 
nal contains signatures of both the black hole QNMs and 
of the backscattering effects. The first ones are triggered 



when the maximum of the rest-mass density crosses the peak 
of the potential at r ~ 3M yielding the largest contribu- 
tion to the waveform at u/M ~ 112. We have verified 
the presence of the black hole QNMs with complex fre- 
quency uji (with i = 1,2, ...) by fitting the waveform from 
this time onward with a superposition of modes of the form 
* (e) = E 2 [ai cos(K(wi)i) + A ninfKCwi)*)] exp[3^)t], 
where the free coefficients ai and Pi are real. The right 
panel of Fig. [2] reports the result of such fit obtained using 
the first five QNMs, and which is essential to properly capture 
the early-time part of the waveform. (Adding higher modes 
does not improve the fit, while considering less modes is not 
enough to reproduce accurately the waveform.) When an- 
alyzed on a logarithmic scale (see bottom part of the right 
panel of Fig. [2]l, it becomes evident that after u/M ~ 130 the 
backscattering effects become dominant and the waveform is 
no longer a simple superposition of exponentially decaying 
modes. 

Additional information on the black hole response comes 
from the study of the energy spectra for shells of different 
initial width a and starting from ro = 15M. This is sum- 
marized in Fig. [3] whose left panel shows the energy spectra 
of the complete waveforms and the characteristic interference 
patterns already introduced in Paper I. Clearly, the interfer- 
ence effects at higher frequencies (i.e., for 2Muo > 0.2) as 
well as the efficiency in the emission of energy via gravita- 
tional waves are increased as the compactness of the shell 
is increased, progressively tending to what is expected for a 
point like particle (cf. the top-right panel in Fig. |4j). Equally 
clear is that the energy spectra do not show any peculiar be- 
havior around the QNMs frequencies (i.e., for 2Mw > 0.35) 
but also that the behavior of the spectra at low frequencies 
does not change significantly with the matter compactness. 
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Figure 3: Left panel: Energy spectra and their typical interference fringes as computed from the waveforms produced by shells inf ailing from 
r\) = 15M. Right panel: The same as in the left panel but only for the ringdown phase, showing the excitation of the 1 = 1 black hole QNMs. 
See text for details. 



As argued in Paper I and in Refs. (HI |23t], interference 
fringes naturally arise in the energy spectrum generated by 
the superposition of two waveforms separated by a certain 
time lag T. Indeed, given two time series A\(t) and ^(t) = 
CAi(t + T), the power spectral density resulting from their 
superposition is given by 

dE 

— oc IAAujM 2 ri+C 2 + 2Ccos(wT)l , (10) 
duj 

where A\(ui) is the Fourier transform of A\(t). As a result, 
such a spectral density will have peaks with a constant spac- 
ing given by Au = 2tt/T. If C = 1, the minima occur at 
dE/dui — and the spectrum is said to have 100% frequency 
modulations. The modulation is always smaller for generic 
values of C. 

The energy spectra displayed in the left panel of Fig. |3]can 
be explained within this general picture. The series of fringes 
is mainly determined by the interference of two bursts of ra- 
diation with different amplitude and separated in time by the 
scattering off the curvature potential at different times. In the 
case of a shell of width a/M = 0.050, the two main pulses 
responsible for this modulation can be distinguished in the 
waveform shown in the left panel of Fig. [2] These correspond 
to the peak produced by the initial motion of and appear- 
ing at very early times (i.e., at u/M ~ 20), and to the one 
emitted when the bulk of accreting matter crosses the peak of 
the Zerilli potential (i.e., around u/M ~ 112). ll48ll The time- 
lag between the two pulses is Au/M ~ 92, which gives a 
separation between the peaks of the spectrum 2MAu> ~ 0.07. 
This value is in good agreement (given the difficulty in unam- 
biguously catch the time when the ringdown starts) with the 
2MAw ~ 0.1 that can be read off from the solid line in the 
left panel of Fig. [3] 



The right panel of Fig. [3] on the other hand, shows the spec- 
tra obtained by performing a Fourier transform of the signal in 
the "ringdown" phase only. The two dotted vertical lines indi- 
cate the n = 1 (fundamental) and n = 2 (first overtone) QNM 
frequencies of the black hole. It is interesting to note that with 
the exception of very wide shells (i.e., with a/M = 1.582) 
which are not able to excite cleanly the QNMs, the maximum 
of the spectra always lies between the two lines, suggesting 
the presence of the two modes in the waveform. 

Another parameter influencing the energy spectra is the ini- 
tial location of the shells rrj. As mentioned above, the largest 
part of the gravitational-wave signal is emitted when the "cen- 
ter" of the shell reaches the peak of the Zerilli potential and 
this will clearly depend on the initial position of the shell. 
Shells that start further away will have longer infalling times, 
larger separations between the first and second peaks in the 
waveforms, and thus smaller separation in the interference 
fringes of the energy spectra. 

In order to reproduce this dependence we have performed 
a number of simulations for shells of fixed initial width 
(a/M = 0.050) accreting from different initial locations rrj. 
Two representative and extreme cases for ro — 10M and 
ro = 70A/ are shown in Fig.|4j whose left panels display the 
waveforms, while the right ones the energy spectra. Overall, 
these plots show a number of interesting and general proper- 
ties of the excitation of black hole oscillations through accret- 
ing matter. Firstly, as the infalling time is increased, the ring- 
down phase is progressively dominated by the non-oscillatory 
tail, essentially as a result of the "spreading" of the shell in- 
duced by the tidal field. This behavior is indeed consistent 
with the analysis carried out in Ref. ll45ll and which pointed 
out that the excitation of the black hole QNMs through a con- 
tinuous flux of infalling particles is made difficult by the pres- 
ence of destructive interferences effects. 
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Figure 4: Infalling shells of dust: waveforms (left panels) and energy spectra (right panels) for two representative values of ro with a/M = 
0.050. For the sake of comparison, the right panels show in dashed lines the 1 = 2 energy spectra for a particle plunging radially onto the 
black hole (rescaled by a convenient factor). 



Secondly, the efficiency in gravitational-wave emission de- 
creases for shells falling from larger distances and is signifi- 
cantly smaller than for point-like particles. This can be easily 
appreciated in the right panels of Fig. [4] when comparing en- 
ergy spectra associated with shells (solid lines) with the corre- 
sponding t = 2 energy spectra of a particle plunging radially 
from the same initial position (dashed line) 04911 . Note that the 
efficiency is in this case of at least a couple of orders of mag- 
nitude larger (the data have been properly rescaled to aid the 
comparison). 

Thirdly, the energy emission is progressively peaked toward 
lower frequencies as ro increases. This behavior can be made 
more quantitative by defining the "characteristic frequency" 



of the energy spectrum as the weighted average 

foj(dE w /duj)du 
Uc = J{dE*>/du)dLj ' (U) 

and by computing how this characteristic frequency changes 
along a sequence of shells of fixed initial compactness and 
falling from increasingly larger radii. A summary of this de- 
pendence is offered in Fig.[5]for shells with a/M = 0.158. 
By following a phenomenological approach, if one defines 
x = log 10 (r /A/) and y = log 10 (2Mw c ) the figure shows 
that the data can be very well fitted by means of a quadratic 
law like y — ax 2 + bx + c with a — —0.461, b = 1.002 and 
c = —0.904. The interpretation of this relation is still unclear 
but certainly deserves further attention. 
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Figure 5: Dependence of the characteristic frequency in the energy 
spectrum on the initial position of the shell. 

As we anticipated at the beginning of this section, an impor- 
tant comment is worth making on the influence of the initial 
amount of gravitational radiation on the waveforms and on 
the energy spectra. Two different approaches are possible in 
this respect and it is important to bear in mind that they do not 
yield identical results. The first one [50] consists in selecting a 
time-symmetric and conformally flat (or non-conformally flat) 
initial profile of the Zerilli-Moncrief function from the solu- 
tion of Hamiltonian constraint. As we discussed in Sec. Ill Bl 
this does not prevent that spurious radiation is produced as 
the evolution starts. The second approach, which has been 
the one adopted here, consists, after the solution of the con- 
straint, in the removal of the spurious burst of gravitational ra- 
diation by evolving the perturbation equations but not the per- 
turbing sources. In this way the initial radiation is allowed to 
leave the computational domain and the evolution can there- 
fore start self-consistently once the initial, nonstationary part 
of the solution has been removed. Figure [6] summarizes the 
impact on the energy spectra and on the waveforms of the two 
approaches for a shell initially at rest at ro = 7.5M, Plotted 
with the dashed line in the main panel is the energy spectrum 
in the case of actual conformally flat initial data involving the 
solution of Hamiltonian constraint only; on the other hand, the 
solid line refers to the spectrum one obtains when the initial 
pulse is allowed to be radiated away. Note the much stronger 
modulation present in the former case and the corresponding 
larger variation in the initial part of the waveform as shown in 
the inset. Note that we would have found larger differences 
for smaller ro and smaller differences for larger ro. Further- 
more, we note that if we start with non-conformally flat initial 
data (that is, throught the solution of Hamiltonian constraint 
with (3 7^ 0) we would obtain larger differences in the early 
part of the waveforms and thus larger modulations in the en- 
ergy spectrum. This behavior is qualitatively similar to what 
discussed in Ref. [23] in the case of particles plunging radially 
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Figure 6: Conformally flat initial data and spurious bursts of radia- 
tion: energy spectra obtained by solving the Hamiltonian constraint 
with /3 = only (dashed line) and when the initial gravitational- 
wave pulse is additionally eliminated (solid line). The shell has 
a/M = 0.158, it falls from ro = 7.5M and the signal is extracted 
at r b s = 500M. The inset shows the corresponding initial part of 
the waveform. 

from finite distance 

We conclude this section by commenting on the amount of 
energy released in gravitational waves as a function of both 
ro and a. This is shown in Fig- [3 which displays in the main 
panel the normalized energy emitted in the lowest multipole 
as a function of the initial location ro and for a shell of initial 
width a/M = 0.158. Note that in the case of a shell this is 
a monotonically decreasing function of the radial distance, an 
opposite behavior to that seen for accreting particles where, 
for conformally flat initial data, the energy has a local min- 
imum and then increases monotonically with ro d], asymp- 
toting the Davis-Ruffini-Press-Price limit In addition, 

the inset in Fig.|7]displays a comparison between the energy 
emitted through the whole waveform (empty circles) and the 
one computed considering the ringdown phase only (filled cir- 
cles). Note that while the black hole ringdown contributes for 
~ 30% — 40% of the total energy in the case of small ro, this 
value goes down to ~ 10% for larger initial distances as a re- 
sult of the progressive loss of compactness experienced by the 
shell during the infall. 



C. Thick accretion disks 

While our model with infalling quadrupolar shells allows us 
to capture some of the essential features of the gravitational- 
wave emission from extended matter sources, such as the ex- 
citation of the black hole QNM ringdown and the presence of 
interference effects in the energy spectra, it nevertheless re- 
mains a useful toy model. 
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Table I: Stable (Do and Di ) and marginally stable (D2) constant angular momentum thick disks orbiting around a Schwarzschild black hole 
of mass A/ = 2.5Mq. From left to right, the columns report the name of the model, the number of radial and polar gridzones used in the 
hydrodynamicsl simulations, the disk-to-hole mass ratio, the polytropic constant « of the isoentropic EOS p — ftp 1 with 7 = 4/3, the value 
of the specific angular momentum I, the position of the cusp r cusp and of the center r ccntCT of the disk, the rest-mass density at the center p c , 
the location of the inner (ri n ) and outer (r out ) disk boundaries, the value of the potential barrier AW, and the orbital period at the center t OI b- 
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Figure 7: Main panel: Total energy radiated in gravitational waves 
(1 = 2 multipole) as a function of the initial location ro of the center 
of a shell with a /M = 0.158. Inset: Total energy versus a for shells 
falling from ro = 15M (open circles). The filled circles refer to the 
energy contribution coming from from the black hole ringdown only. 

To improve on this and to examine a scenario which 
is astrophysically more realistic, we will now consider the 
gravitational-wave emission resulting from the accretion of 
geometrically thick disks (or tori) orbiting around a non- 
rotating black hole. Such systems are believed to form in 
a variety of different ways, such as during the last stages of 
gravitational collapse or in the merger of binary neutron stars. 
In addition, if compact enough and undergoing oscillations, 
these systems generate gravitational-wave signals within the 
sensitivity curve of ground-based interferometers lfl8ll26ll27ll . 
Hereafter, and as in the case of dust shells, we will assume that 
the mass of the torus /i is much smaller than that of the black 
hole. 

Detailed descriptions on how to build equilibrium configu- 
rations of barotropic thick disks orbiting black holes are given 
in Refs. f24l l25ll . but we here recall that these objects have 
traditionally been described as obeying a polytropic equation 
of state p = Kp 1 , with 7 = 4/3. Because in pure orbital 
motion, the fluid four-velocity is given by u a — (u f ,0, 0, vr") 



and it describes a non-Keplerian rotation around the central 
black hole, with a specific angular momentum £ = —u^ju t 
distribution which is essentially unknown. 

Table U lists the main features of the tori considered in our 
simulations, such as the values of the specific angular momen- 
tum I (assumed constant throughout the disk), of the potential 
gap AW between the inner edge of the disk and the cusp, and 
of the mass of the disk fi. In order to induce a non-trivial dy- 
namics in these otherwise stationary disks, the initial models 
have to be perturbed in some way, such as by adding a ra- 
dial velocity v r = 77t^ ph , where 77 is a parameter and w^ ph 
is the radial velocity of the spherical stationary atmosphere 
surrounding the torus 12611 . The rest-mass density of the at- 
mosphere is low enough not to influence the dynamics of the 
tori (its density is several orders of magnitude smaller than the 
maximum density at the center of the disk). 

As in the case of infalling dust shells, the initial data for the 
perturbation fields ^% is obtained either through the so- 
lution of the Hamiltonian constraint (for even-parity perturba- 
tions) with conformally flat initial data, or assuming stationar- 
ity, i.e., through Eqs. (O with df^'fj ^ = 0. Clearly, in either 
case this initial data is not consistent with the hydrodynami- 
cal sources once they are perturbed. This mismatch inevitably 
introduces a unphysical initial burst of gravitational radiation 
that is however easy to distinguish and remove from the anal- 
ysis. 

We note that the accretion of matter onto the black hole 
can occur on different timescales depending on the perturba- 
tion and we will consider here the two limiting situations that 
could be encountered in astrophysical scenarios. In the first 
case we will consider a perturbation which is large enough 
to cause a runaway accretion of the disk onto the black hole; 
we will refer to this as the "hypercritical accretion" scenario 
(see Sect. IIVC II ). In the second case, on the other hand, the 
perturbation is smaller and it will induce a series of quasi- 
periodic oscillations, each accompanied by an episode of ac- 
cretion onto the black hole; we will refer to this as to the 
"quasi-period accretion" scenario (see Sect. lIV C~2l . 



1. Hypercritical Accretion 

As mentioned above, simulations of hyper- accreting disks 
can be performed by simply choosing a sufficiently large ini- 
tial value of the radial velocity perturbation. Under this con- 
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Figure 8: Model Dq. Left panel: time evolution of the equatorial rest mass density profiles for a velocity perturbation r\ = 0.3. Right panel: 
time evolution of the center of mass in the equatorial plane. The disk disappears in the black hole at t » 0.75 i or b. 



dition, the centrifugal barrier cannot counteract the injection 
of additional kinetic energy and the whole disk is "pushed" 
towards the black hole, rapidly accreting onto it in less than 
one orbital period. 

Hereafter we will concentrate on the dynamics of model 
Dq, to which a velocity perturbation with r) = 0.3 is added. 
For this model the hydrodynamical grid runs from r m ; n = 
2.03 M to r max = 16 M and it is covered by 300 x 150 
grid-points, geometrically spaced in the coordinate r, with 
Ar m i n = 3 x 10~ 4 and Ar max = 0.2. The range of the 
r* ID grid is r* G [-2200 M, 4500 M], to guarantee that the 
outer boundary does not influence the slope of the late time 
power-law tail, and is covered by~ 3 x 10 4 cells. 

The left panel of Fig.[8]shows, in terms of the tortoise coor- 
dinate, five snapshots of the evolution of the rest-mass density 
on the equatorial plane. Note that as time proceeds, the disk 
is compressed, the density increases and matter quickly starts 
falling onto the black hole. This continues until the entire disk 
has been accreted. A summary of this is shown in the space- 
time diagram in the right panel of Fig. [HJ which presents the 
motion of the projection on the equatorial plane of an effective 
"center-of-mass" of the disk, defined as 



r pr 2 dr 



(12) 



where g rr = (1 — 2M/r)~ 1 (for convenience the spacetime 
diagram is shown using the tortoise radial coordinate r*). Cor- 
respondingly, the inset in Fig. [8] depicts the time evolution of 
the mass accretion rate for model Do. 

The vertical dashed line in the right panel of Fig. [8] is ap- 
proximately the location of the peak of the Zerilli potential, 

r pcak _ 19QM (r pcak _ 3 j M) and can be used tQ Men . 

tify t ~ 0.44 i or b as the time at which the center of mass 
crosses the peak of the potential. Furthermore, this panel 
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Figure 9: Model Do: gravitational waveform (extracted at r b s = 
200 M). The presence of two ringings is evident. The first one is 
related to the gravitational-wave content present in the initial profile 
of the Zerilli-Moncrief equation. The second ringing (magnified in 
the inset) is determined by the object crossing the peak of the poten- 
tial. The dashed line, proportional to u~ 7 , shows that the late-time 
behavior of the waveform is properly captured. 



also shows that the center of mass leaves the hydrodynami- 
cal grid (r™ in ~ —6.35 M) at t ~ 0.75 t olh and from this 
time onward, the gravitational-wave signal is dominated by 
the QNMs of the black hole. 

The gauge-invariant gravitational waveform produced by 
the hypercritical accretion has a dominant I = 2 character lplll 
and is displayed in a logarithmic scale in Fig. [9] as measured 
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Figure 10: Model Di: time evolution of the projection of the center 
of mass in the equatorial plane. After an initial transient, the dynam- 
ics is characterized by constant amplitude quasi-periodic oscillations. 
The inset shows the mass accretion rate. 

by an observer located at r Q b s = 200 M. Besides the ini- 
tial burst of spurious radiation related to the initial data, and 
which is no longer present by u ~ 0.2i or b, the waveforms 
exhibits two main features. The first one is determined by 
the infalling matter as its center of mass approaches the max- 
imum of the potential at u ~ 0.4t or b- The second one, in- 
stead, starts at u ~ 0.7i O rb and represents the ringing re- 
sulting from the perturbation experienced by the black hole 
through the rapid accretion of the torus (see the mass accre- 
tion rate inset in the right panel of Fig. [8]). This becomes ap- 
parent when comparing the features of the waveform in Fig. [9] 
with the matter dynamics in Fig. [8] The ringdown phase ends 
with the usual power-law tail; this is highlighted by the dot- 
dashed line in Fig. [9] which is proportional to u~ 7 and shows 
that for u > 3.6£ or b the waveform has only a quadrupolar 
nature decaying as oc u~( 2e+3 \ 

In summary, the analysis of the process shows that a very 
rapid and hypercritical accretion of a compact distribution of 
matter onto a black hole produces gravitational-wave emis- 
sion that is qualitatively similar to that of either relatively 
narrow dust shells plunging from a large distance or wide 
dust shells accreting from a small distance. Furthermore, the 
gravitational-wave signal shows oscillations that reflect the 
excitation of the black hole QNMs (especially the fundamen- 
tal one) but these oscillations are in general so weak that the 
non-oscillatory tail determined by the long-range properties 
of the scattering potential soon dominates the signal. 



2. Quasi-periodic accretion 

A phenomenology complementary to the one discussed in 
the previous section for a hypercritical accretion can be stud- 



ied when the accretion takes place in a quasi-period fashion. 
Indeed, there are at least two different ways of triggering a 
quasi-periodic accretion and we will discuss them separately. 

The first possibility has already been discussed in 
Sect. |IV Cl and consists in adding a perturbation in the radial 
velocity of the torus parameterized in terms of an analytic ac- 
cretion solution. As a concrete example we will now concen- 
trate on the dynamics of model D\, describing a disk with a 
comparatively high value of the specific angular momentum, 
I = 3.80, and subject to an initial radial velocity perturbation 
with jj = 0.2. 

Figure [10] shows the motion of the center of mass on the 
equatorial plane, indicating that after an initial transient, a 
phase of quasi-periodic, (almost) constant amplitude oscilla- 
tions follows, in which the disk periodically approaches the 
black hole before the centrifugal barrier pushes it back, past its 
original position. As the torus nears the black hole, part of its 
matter is spilled through the cusp, resulting in a quasi-periodic 
accretion of matter onto the black hole; this is shown by the 
inset of Fig. [TOl which reports the evolution of the mass accre- 
tion rate (see Ref. 112611 for a more detailed discussion). Note 
that the accretion is essentially shut-off as the torus moves 
away from the black hole. 

This quasi-periodic dynamics of the torus is clearly im- 
printed onto the gravitational-wave signal and this is shown 
in the left panel of Fig. QT| where the I = 2 gauge-invariant 
waveform (extracted at r Q b s = 500 M) is superposed to the 
signal extracted using the quadrupole formula (SQFi). Note 
that the very early time part of the gauge-invariant signal has 
been removed to avoid influences coming from the burst re- 
lated to the initial data. The waveforms show a small burst 
at u w f or b, that corresponds to the initial accretion of mat- 
ter, followed by regular oscillations that mirror the motion of 
the torus in the potential well. Some differences between the 
gauge-invariant and the SQFi waveforms are recognizable in 
the initial stages, with the first one showing highly damped 
high-frequency oscillations for u w 2 f or b superposed to the 
main, quasi-periodic signal. Since this high-frequency pat- 
tern is absent in the quadrupole waveform, we conclude that 
it originates from curvature backscattering of gravitational 
waves emitted in the initial infalling phase. 

We also note that the differences between the waveforms 
computed through the gauge-invariant formalism and the 
quadrupole formula using the quadrupole tensor defined by 
Eq. (O become larger if the initial velocity perturbation is in- 
creased. This is simply due to the fact that the SQFi is valid 
in the weak-field and slow-motion approximation and thus 
ceases to be accurate when the dynamics of the oscillations 
is nonlinear, with the velocities becoming relativistic. On the 
other hand, the velocity-dependent terms in the source func- 
tion of the Zerilli-Moncrief equation are obviously accurate 
also in a relativistic regime. 

The energy spectrum for this signal is shown as a solid 
line in the right panel of Fig. [TT] and refers to the signal from 
u/t 



orb 



> 



1 in order to eliminate the unphysical initial burst. 
The spectrum shows the presence of a broad-band peak at 
low frequencies, related to the highly damped oscillation in 
the early time part of the I = 2 signal, together with nar- 



13 




Figure 11: Model D\, Left panel: gravitational waveforms (extracted at r b s = 500 M) versus retarded time. h + polarization (even-parity) 
superposed to the waveform extracted using the SQFi (dashed line). Right panel: Energy spectrum for the I = 2 and I = 3 (odd-parity) 
modes. The highly damped oscillation in the early part of the I — 2 waveform (from $20 ) resu lt s m a broad-band peak which is not related to 
the QNMs frequencies of the black hole. 
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Figure 12: Model D2. Left panel: the motion of the center of mass on the equatorial plane and the mass accretion rate (in the inset). Right 
panel : the qudrupole gravitational waveforms. The gauge-invariant (extracted at r Q b s = 500M) and the SQFi waveforms are shown on the 
same plot for comparison. 



row peaks in the harmonic ratio 1:2:3:..., which in- 
stead mirror the non-linear coupling of modes in an oscillat- 
ing torus (see Refs. 11281 12911 for a detailed discussion of the 
eigenfrequencies and their astrophysical implications). Such 
peaks are present in both the I = 2 and I = 3 energy spectra 
(dashed line in Fig.fTTTi. The latter is generated by a waveform 
rh 3 ® / (cos 9 sin 2 9) whose amplitude is roughly two orders of 
magnitude smaller than the I = 2 one. It is interesting to 
note that while for the I = 2 multipole the first frequency of 



the disk is the dominant one, the second frequency has the 
largest power for the I = 3 multipole. The energy spectrum 
shows no distinctive signature at 2Mlo ~ 0.7473 that could 
be related with the excitation of the QNMs of the black hole 
(cf., right panel of Fig. [3j. This is in contrast with the results 
of Ref. 11911 . which have simulated a similar source of per- 
turbations for the black hole and have indeed found a minute 
high-frequency contribution in the energy spectra related to 
the emission from the black hole ll52ll . It is presently unclear 
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what is the origin of this difference, but it is likely that the 
hybrid approach proposed in Ref. 01911 . and which combines 
a solution in the time-domain for the hydrodynamics and one 
in the frequency-domain for the perturbative equations, is bet- 
ter suited to extract the extremely small contributions coming 
from the black hole. 

A second possibility of triggering quasi-periodic oscilla- 
tions in the disk, and hence a quasi-periodic accretion, is by 
an instantaneous reduction of the specific angular momen- 
tum. The main difference with the dynamics discussed before 
is that this reduction produces a continuous mass accretion, 
which is however modulated quasi-periodically as the torus 
approaches or moves away from the black hole. This behav- 
ior is shown in Fig. [12] for model D 2 , in which the specific 
angular momentum is reduced from I = 3.80 to I = 3.72. 
More specifically, the left panel of Fig. Q~2] shows the motion 
of the center of mass on the equatorial plane and the evolution 
of the mass-accretion rate in the inset. Because of the intrinsic 
reduction of the centrifugal support and the resulting contin- 
uous accretion, the oscillations are quasi-periodic in time but 
with decaying amplitude, so that by t/t OI b ~ 50 they have es- 
sentially disappeared and the torus accretes at an almost con- 
stant rate, with the center of mass progressively approaching 
the black hole. 

This different quasi-periodic dynamics of the torus is also 
imprinted onto the quadrupolar waveform and this is shown 
in the right panel of Fig. Q~2] which presents both the gauge- 
invariant signal as measured by an observer located at r ^ s — 
500 M (solid line) and the one computed with the SQFi 
(dashed line). Note that the differences between the sig- 
nals are very small initially, but they grow with time and 
by t > 15£ or b, when the angular momentum loss induces a 
global radial motion of the disk, the mean value of the gauge- 
invariant signal progressively drifts away from zero. This 
different behavior can be understood simply if one bears in 
mind that during this stage the disk is undergoing an almost 
steady-accretion phase. As a result, while the disks mass- 
quadrupole is nonzero and large, its time variations, and hence 
the corresponding gravitational-wave signal as computed in 
the SQFi, are essentially zero. The perturbative approach, on 
the other hand, is able to capture the perturbations induced by 
this steady-state accretion and this is reflected in the gradual 
and secular increase of the gravitational-wave signal shown in 
the right panel of Fig. [12] Of course, this behavior is the same 
already encountered in the case of accretion of dust shells and 
it is, once again, generated by curvature backscattering off the 
"tail" of the potential. 

Note that during the final stages of the accretion, i.e., for 
t > 60 t or },, the motion of the center of mass exhibits an expo- 
nentially rapid motion toward the black hole which, however, 
is not reflected neither in the mass-accretion rate, nor in the 
gravitational-wave signal. The reason for this is that, in prac- 
tice, the rest-mass contained in the torus in these final stages 
is extremely small and the rest-mass density is so diluted that 
it becomes comparable to that of the atmosphere surrounding 
the torus. As a result, the effective perturbation induced onto 
the black hole in the final stages of the accretion is vanishingly 
small. 



In summary, the analysis of this process shows that the 
gravitational-wave signal produced by a quasi-periodic accre- 
tion can be rather different, depending on whether the ampli- 
tude of the oscillations is (almost) constant in time and the 
accretion rate is also periodic, or the oscillations have de- 
creasing amplitude and the accretion rate is constant with a 
periodic modulation. In the first case, the gravitational-wave 
signal will reflect faithfully the dynamics of the matter with 
constant-amplitude waveforms averaging to zero and a spec- 
trum showing peaks at the eigenfrequencies of the oscillating 
matter. In the second case, on the other hand, the waveforms 
will have a decreasing amplitude and will not average to zero 
as a result of the underlying continuous accretion. 



V. CONCLUSIONS 

By performing numerical simulations that combine the so- 
lution of the nonlinear hydrodynamics equations with that 
of the linear inhomogeneous Zerilli-Moncrief and Regge- 
Wheeler equations, we have studied the features of the 
gravitational-wave signals generated by the accretion of mat- 
ter onto a Schwarzschild black hole. 

As extended and accreting matter-sources we have consid- 
ered quadrupolar shells of dust falling radially from a finite 
distance, as well as geometrically thick disks undergoing ei- 
ther bursts of hypercritical accretion or quasi-periodic oscil- 
lations. In both cases we find that the gravitational-wave sig- 
nal is not a simple superposition of the black hole QNMs and 
that the latter cannot be found in the energy spectra at times. 
Rather, we find that quite generically the signal contains im- 
portant contributions coming from radiation scattering off the 
tail of the curvature potential and producing a characteristic 
pattern of interference fringes in the energy spectra. While 
the relevance of this contribution differs according to the spe- 
cific source considered, it is generically present as long as the 
source of perturbations is extended and the scattering potential 
does not have an exponential decay with radius. This conclu- 
sion, which was already reported in previous studies involving 
simpler sources, has been here confirmed unambiguously by 
studying the scattering off a fictitious potential, the Poschl- 
Teller potential, which however decays exponentially with ra- 
dius. 

These generic properties of the gravitational-wave emis- 
sion coming from black holes perturbed by extended sources 
represent important differences with respect to correspond- 
ing properties of signals produced by very compact sources, 
such as point-like particles. Of course this conclusion pre- 
vents from the derivation of a simple and generic description 
of the gravitational-wave signal which would be independent 
of the properties of the perturbing source, but it also opens the 
exciting perspective of deducing many of the physical proper- 
ties of the source through a careful analysis of the waveforms 
produced. 

Overall, the results presented here make us confident that 
the black hole QNM contributions to the full gravitational- 
wave signal should be extremely small in generic astrophysi- 
cal scenarios involving the accretion of extended distributions 
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of matter. On the other hand, it should also be stressed that the 
time-domain analysis carried out here may not be the most ac- 
curate to extract the contributions coming from the perturbed 
black hole when these are several orders of magnitude smaller 
than those coming from the source itself or from the backscat- 
tering off the potential. In these cases, however, a hybrid ap- 
proach such as the one proposed in Ref. lfl9tl . combining the 
solution in the time-domain for the hydrodynamics with one 
in the frequency-domain for the perturbative equations, may 
be better suited. 
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